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Some Application of Modern Mathematical Statistic to the 
Lunar Tidal Analysis of Ionospheric Height. 


By K. MAEDA and M. IZAKI 


(Electrical Communication Laboratory, Ministry 
of Telecommunications.) 


Summary 


The paper deals with some problems of lunnar tidal anlyses of 
ionospheric height, which are treated from the standpoint of mathe- 
matical statistics. The semi-amplitude of lunar semi-diurnal tide of 
ZF: (height of maximum electron density of Fs region) analyzed 
from Japanese data is tested by statistical hypothesis and its results 
are described. Next the error caused by a rough scaling of height is 
estimated in connection to scaler’s skill (or carefulness), and the results 
are given. 

Some statistical measures are introduced for the criterion of 
the validity of analyzed results in the above treatments, that is, the 
level of significrnce for the former case (test of statistical hypothesis 
by non-parametric method) and the-coefficient of confidence for the 
later case Joudess of error). 

I. Introduction 

(with particular reference to the work by Chapman and Bartels) _ 

The movements in the ionosphere, such as wind, circulation, tidal motions etc., 
have recently been considered as one of the most importants problems of the iono- 
sphere physics. The lunar tidal phenomena of the ionospheric height have been 
revealed by many scientists on the basis of various analyses of observed results and 
it is reported that amplitude of the semi-diurnal lunar tide of the ionospheric height 
will be a few kilometers. 

In the harmonic analyses of ionospheric data careful considerations should be 
paid to the observed materials themselves, their selection and arrangement for analyses 
and the method of analysis, and fairly elaborate works are necessary to evaluate the 
‘errors and to test the periodicities. Regarding these problems we can find and 
extensive study in “Geomagnetism” vol. II. by S. Chapman and J. Bartels. Although 
this book covers every item of importance for the harmonic analysis, it shall not be 
“unnecessary to treat a few problems of analyses on the basis of modern mathematico- 
statistical conception, the methods oe which will se ety described in mentr est to 


those written in Geomagnetism, — 


In this book the investigation of periodicity is discussed (§30 in p. 593 and 
equation (86) in § 25), but it is necessary to know a parameter in population, that is 
the variance, which is put approximately equal to the sample variance in the case of 
a great number of data. The present authors made the test of statistical hypothesis 
by using the non-parametric method, such as F-distribution. 

Also in the book (p. 581 and others) the effect of observation errors is discussed, 
and we will discuss the effect of scale in connection to scaling skill (or carefulness), 
which is different from that described in the book. 

Before making the harmonic analysis, or applying the test of statistical hypothe- 
sis on periodicities and estimating the errors, it is necessary to make the process 
stationary (for example, to eliminate trends) and to investigate the elements, which 
construct the process in population (to make the correlogram), and finally to perform 
the test of statistical hypothesis on these investigations. The quasi-persistence intro- 
duced in the book corresponds to the trends and several persistence curves (Fig. 12 
p. 589) involve the same sense as the correlogram. But these investigations will 
become complete after the test of statistical hypothesis on. themselves. 

The problems of stationary process and correlogram are not described here in 
detail, but we want to point out here some problems from the consideration of the 
current observational practice of the ionosphere and to discuss some of them statis- 
tically in the following sections. 

Problem 1. Even if we can measure the ionospheric dead in a sufficient 
accuracy (i. e. height mark of sufficiently small step) and on lunar time basis (i. e. lunar 
hourly or semi-hourly base), it is necessary to test the significance of the analyzed 
result statistically, as sampling errors are involved. 

Problem 2. The scaling practice of the ionospheric height of current use in 
the world is 10 Km step. Such quantization of the height record is comparatively 
great and is seriously connected to the validity of the analyzed results. In this case 
the skill (or carefulness) of scalers should be taken into account. 

The problem is restated as follows. How long should be the period for the 
analysis (a month, a season or a year), or how much should be the number of observed 
materials for the analysis, in order to obtain sufficiently small errors in conriection 
to the quantization of scaling skill (or carefulness) ? 


Problem 3. In most cases of the lunar tidal analyses we have to usé asta 


observed on solar time basis and the difference of time basis (solar and gee will 
exert errors to the analyzed results. EN 

In the first problem” we can apply the test of statistical hal by non- 
parametric method and obtain a definite level of significance for each individual resuit 
of analysis. In rhe second problem we can determine an average error with a certain 
definite coefficient of confidence for each individual results siete from a certain 
number of observed values (sample size). — noida 1 ot inersosings 


In the following the writers describe some typical ittdvinn of sis above two- 
pe Dieca which are treated in the light of mathematical statistics, The third problem — 


is left for the future study. 


If. Errors and distributions of observed values. 


Jonospheric data of Japan were analyzed by Matsushita and it was reported by 
him“ that the semi-amplitude of a half-day period (lunar-time) of .Z,.F2 (height of 
maximum electron density of Fz region) was about 3 Km. In his analysis the errors 
caused by fixing a reading scale of the ionospheric height at 10 Km step and by the 
difference between lunar and solar time bases are involved. Even in cases of accurate 
measurements performed with a sufficiently precise reading and on lunar hourly base, 
sampling errors will be involved so far as the materials for the analyses are observed 
data. Therefore the test of statistical hypothesis or the statistical estimation must 
be carried out in order to account the obtained value (such as 3 Km) as the real 
effect of tidal movement. 

Before the harmonic analysis on observed data is performed, the several steps 
which are described in the preceding section (to make the process stationary and to 
investigate the process by making correlogram etc.) must be taken, but the descrip- 
tion on these points is omitted here, and we assume that the necessary data w#(¢) for 
harmonic analysis is the sum of errors é(f) due to height scaling practice and x(#) 
which should correspond to the phenomena represented by the equation as follows. 


z= Ka, cosd,£+6,sin4,t)+v(4) = de, cos (4,é+¢,)+0()=m(t)+v(#) (1) 


v(t) is a random element and there will be no objection in taking it as Gauss- 
ian distribution N(o0,o°) in our case. Also we assume here that the material can be 
arranged’so as every solar cyclic term is eliminated. 

Thus, we obtain as the equation of abserved value z (#), 


u@)=x (+d H=m()+vO+0O 


and this will show a discrete distribution with the minimum interval of the - height 
mark (say 10 Km). In any case the observed value is more or less discrete. The 
discreteness will not be the problem, but the point is the relative comparison between 
the minimum unit of height scale and the state of distribution (here, 0”). In our case, 
it seems that the scale is fairly rough to take w(¢) as continuous, judging from the 
distribution of u(é). Even if we can take as continuous, it deviates more or i io 
from the Gaussian distribution owing to the errors (06) of scaling. 

... But considering from the practical aspect, it will be better that we first put 
u(t) as N (m (i), 0?) and estimate the level of significance in the test of statistical. 
hypothesis as discribed in the following section, and after that we take the actual. 
level as a little higher than the level based on normol distribution. - 


a Significance, of Harmonic. Analysis. igo peer yt) ere er ee 
; We put the ‘observed value # (4) as N(m (8), o*) tentatively, sanstleting’ the odie 
cussions : in the. preceding section. elt © we Beg eg erat ee ae 


Then the equation 


we 
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indicates the F-distribution with the degrees of freedom (2, M—2h—1) based on the 
normal. regression theory*. In equation (2) N is a number of observed data, and 3} 
working on m% and uw? is the sum of all data, and >} of the last term in the denomi- 
nator is the sum covering every semi-amplitude calculated for a definite number (/) 
of lunar harmonics that were supposed to exist a priori. By giving suitable values 
to a, and b, in the numerator of equation (2), we can perform a test of statistical 
hypothesis on the semi-amplitude under consideration. It is worthy of notice that in 


this equation (2) no parameters in population are involved. 


Here we employ the method of null-hypothesis, that is, we calculate the equa- 


tion (2) by putting c,=0 (a,=0, b,=0).. The value of F of equation (2) for this case 
can be calculated from the data and we put it as /. 

On the other hand from the table already made, we know alemees of F*y—21-1 (a) 
with the degrees of freedom 2 and N—2h—1 and for the level of significance a. And 


in the case when 
Fy > F’ y_2n—1 (@) (33 
we must reject the hypothesis c,=0 with the significance level @. Therefore it can 
be mentioned that: c, is significant, but this mention may not be right with = rate 
of a. ee ye 
In the following we apply the test of statistical hypothesis on some of the results 
of half-day period (lunar time) in the harmonic analysis worked out by Matsushita. 


pi following table I shows three examples reported by Matsushita. 
Table I 


In order to perform the 


_semi-amp. (Km) 
Period 


(A) Feb. 10-Mar. 9, 1947 


3.0 1.4: 
(B) Mar. 22-Apr. 19, 1947 3.0 2.4 
(C) 369 days. from Jan. 6, 1947. 2.7 0.1 


Toh Wich il 


test, it is necessary to know all 
the actual data unmanaged after 
the observation and calculate all 
the semi- amplitudes. But we 


have only the mean values in 
- each lunar hour throughout one 


month (29 days) or one year (369 days) at hour hand. Therefore we treated the 


problem as follows for a preliminary trial. 


We put N of (N—2h—1) in the equation (2) as 2429 in the case of oie month. 


and as 24x 369, in the case of one year and calculated v= ue = fos Ue y by means 


* See for example, S.S. Wilks: Mathematical statistics. 


1943. 


Chapter 8, Princeton Unive ner eee z. 


of ee ee eee a ee 
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of 24 data (here N=24) and later multiplied it by 29 in the former case and by 369 
in the later case in order to obtain. 
The values of Fy in the cases (A), (B) and (C) and F2y~2n—1 (a) are given in the 
following table IL. 
Table II In the case (A) 3 Km as 
- 2 (%) Sm the semi-amplitude of lunar semi- 
diurnal tide is significant with 
the level of significance of 10% 
an ; : and in the case (B) not signifi- 
(B) 1.4 6.91 4.64 | 3.01 | 2.30 | 1.61 cant even with 20%. In the case 


01 | 1.0 | 5.0 {10.0 | 20,0 


_ (A) 2.9 6.91 4.64 3.01 2.30 1.61 


(Cc) | 11.35} 6.91 | 4.60 | 2.909 | 2.30 | 1.6, (©) it is significant with the level 
ieee | of 0.1% 

Judging from the considerations discussed in the preceding section and the data 
used, it is necessary to estimate the level of significance as a little higher than those 
listed in the table IJ. In the case (B) the level-of significance will be heigher than 
20% and in the case (A) the level of 10% may not be safe. Therefore it will be safe 
to consider the level of significance as more than 20% for both cases. But F% in the 
case (C) is sufficiently larger than 6.91 corresponding to the level of 0.1%, and it can 
be mentioned that the results of case (C) (2.7 Km) is significant with at least 1 % of 
the significance level. 


A doubt remains still on the significance of 2.7 Km’ as the semi- -amplitude of 
lunar semi-diurnal tide, considering that it was an analyzed result by using the data 
observed with a fairly great scaling step of 10 Km. In the following section we 
investigate the problem of errors caused oy scaling. se 


IV. Influence on a harmonic analysis reed by the errors attributable to scaling. 

1. Errors of a harmonic analysis 

When we carry out a harmonic analysis, such as a periodgram analysis, by 
using 2 (2), a periodgram h*y is given as follows, 


om NB (t)48 (O)} cos at J +[ SoH (+6 ()}sin 2 i} | 
= (2 0082 +(3} x sin A ¢)? 422 (4) 


The first term in the above equation (4) is a periodgram in case of lacking error (0), 
and this is the actual value we need. And Z? is the error of a periodgram which is 
attributable to 6. If we can mention with a large probability that Le is small enough 
in comparison with ¢, in eq. (1), it must be recognized that the analysis using #(#) is 
not senseless. Although we must know distribution of Z accurately for the estimation 
of Z, we try to estimate Z by calculating only E(Z’) or E(Z') (E means average). 

After comparing the various methods which can n be employed the following one 
will be the most suitable for the present. problem.) = 


4 Sie 3 re See, CAE ae 2 i ” Hs Sand BY 


6 
The equation (5) shows the relationship between the minimum of the estimation of 
error (Z) and the given value of the lower limit of the confidence coefficient (<). We 


neglect the equations corresponding to £(Z') or the. similar ones. 


The calculated results of E(Z”) is 
BZ) =A SUE OIE CO x(09) igh 


a 


2. Scaling and scaling skill (or carefulness) 
We have to investigate the problem of scaling before we proceed forward. 


Though the minimum step of height scale is 10 Km in our case, it is generally put 
as k (Km). 
i) The most skilful Yi le will record the height w(t) as kp, when 


(p: integer) 
OE a 
aa 3 pStshbts 2 (x: actual height) 


ii) In the case of very unskilful observer the probability with which he records — 
the height a(t) as kp will be equal to the epeecthon ge gy SH he records u wh poe as 
kp+k when kpsaxshp th : 5 een 

That is 


a Fe 


pr (ae kp)= pr(u= p+ b= 


sn 


ii) But it is supnosee that the Acttal cases will be intermediate between the 
: above two situations. a pee pet S es PRE 
This case is represented by : aod icin: Ss a 
r. etuo 2b fiserwif-deres ten’ - 
Pr (u=Kb)=6 (0, - EI 
cada — 


wicks 7) (2) is is an Sonus function of x, and we have to ae PE _sever, 
zi conditions on ¢ (x), considering the actual scaling practices. In general g(a) cont 

oe few parameters which must be determined from the actual scaling. In the a 

; a) Sas  @ and Gil) above men Pade 

Sh Phe a the special cases Mont 


moles ct feavi meter, and oes when 9( 
Z ii 


RE Om a 1 


of a random element v (Z). 
Next we obtain the other equation ei 


E(dx)= 20°S}(— res od =e 


where & is the unit of scale and o indicates the order of a variation. So it may. be 
possible to assume as k<o, and we can obtain the following results with a relative 
error smaller than 10-, 


eae sie 
EO)= 7, oval | (7) 
E(ox)=0 ~ eid 


When the = aad of g(x) isa triangle, we can obtain the results with errors of 
the same order _as above, 


to 


k 


en (9) 
0 


ie 


E(&)= 
a (6x)= 


: 


m2 ‘Estimation ‘oe 
After substituting the pairs of equation (M, (8) and (9), (10) respectively. in ‘the 


equation (6), we get E(Z) and then we calculate the es (5) by using | newly 
obtained Ez yy, The results are as follows, 


In the case (i) : gurrst : ae a 
ne Se (11) 

2 { ean" ; a 

In the case, when the sharp of. ne ig tnttianglé) sesnpe et orb daass os 
somal oil to apssin: orrr <3 bas a1 * i ae Siecrcint Bas ee sens oi (a2 
ee [Tee : 3 a Per sedtiennie) dressed Nl 


at By the eeation (11) and 1 (2), we calcurate the upper limit of os corresponding 


to the lower limit of a confidence coe en = e, putting k=10 Km, and get the results 
as shown in the following table TM 2 Hloignt..L bas. asnrqsH? 2 mes 


ri; JAS “gitioagr “idinom: srltpapie et 1sx0- : blest. M bas bon a co 


oxailg 


abe 
pS SR- 


M2 


ots Pade seca. 


aus eed (= 10Km) Satie tbissesH 


8 
will be fairly large even for the comparatively low coefficient of confidence as 80%. 

V. Conclusions i i . 

We: discussed about the significance and the error concerning the results of 
lunar tidal analyses of ionospheric height and described at first the test of statistical 
hypothesis as applied ‘to some examples of harmonic analyses reported by Matsushita, 
assuming that the observed value z(t) has Gaussian distribution N(m/(2), 0°): Although 
u(t) should be taken to deviate from this distribution owing to the scaling error, it 
is preferable to estimate the actual level of significance as larger to Some extent than 
the calculated level based on normal distribution. Three examples of test for period 
of one month and one year for the lunar tidal analyses of ZmF, are presented. 

Then we investigated errors, which are to be attributed to the scaling, that is, 
minimum step of (k=10 Km) height mark and scalers skill (or carefulness) and pre- 
sented some results of calculated errors. — 

The estimation of error (Z) is proportional to the width of scale (k) and inversely 
proportional to the square root of sample size (V) under some assumptions. In this 
connection the method of expressing the degree of skilfulness (or es 
was discussed. ae 2 ee me : 

The level of significance i in n the test of Paes) rapes and A errors with 
a coefficient of confidence will serve as definite quantitative measures s for criticizing 
the validity of the pie a 2 
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On the Charge Distribution in Voleanic Smoke 


By G. ISHIKAWA,* M. KADENA* and M. MISAKITt 


*Meteorological Research Institute, Tokyo and }Magnetic Observatory, Kakioka 


Abstract 


We made a simultaneous observation of the electric potential 
gradient and the space charge under the volcanic smoke of the new 
crater of Volcano Azuma in Fukushima-prefecture from August 13th 
to 19th 1950, to investigate the charge distribution in the volcanic 
smoke. The remarkable points which we obtained by the analysis of 
this observation are: 
(i) The charge of the smoke of Volcano Azuma is positive, and 
the density of which is estimated at about 5 e.s.u./m’ near the center 
of the smoke. 
(ii) A thin layer of the negative charge is induced around the 
cylinder shaped smoke. The density of the induced charge is about 
1/10 of that of the positive charge near the center of the smoke. 
1. Introduction 
Many investigations confirm the fact-that the volcanic smoke has electric charge 
[1], [2], [3], [4], [5]. Most of these investigations, however, are conducted on the basis 
only of the electric potential gradient at the earth’s surface under the smoke, so that 
it is inevitable that an arbitrariness should enter into the deduction of charge dis- 


tribution in the smoke. For example, even if the variation of the electric field under 
the smoke is negative, we cannot say 


that there is no region of positive 
charge in the smoke. | 
To confirm the charge distribu- 
tion more in detail we must carry out 
the simultaneous observation of the 
potential gradient and the space charge, 
which we made actually at Volcano 
Azuma. Volcano Azuma had been dor- 
mant until the recent eruption on Feb. 
10th, 1950, since the last eruption of 
1896. The recent eruption occurred 
without any symptom, followed by a 
~ second smaller one on the 19th of the 


10 


same month. It was observed that large quantities of ashes, pieces of lava and rocks 
were scattered as far as 600 meters from the crater. When we made the observation 
there, Volcano Azuma was already in a quiet condition as before. Its smoke was 
white, looked like a steamjet, and was supposed to contain a very small amount of 
solid particles. The state in such smoke must be greatly different from that in black 
smoke. The observing point O was in the leeward of the crater and the smoke flowed 
over it almost continuously, and frequently it came down to the surface of the slope, 


covering the observing instruments. The smoke remained in general the circular 


A 


4058 ~(0.5.tt,/n) 
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cylinder shape, whose diameter was estimated to be about 30 to 50 meters, except 
when it swept down to the ground. ° 


2. Observing Instruments and the Results Obtained 


The instrument used in measuring the electric field is the portable electrometer 
with a radioactive collector mounted on a tripod. Owing to the rocky, complex con- 
figuration of the ground, it is difficult to make plane reduction for the observed values 
of the potential gradient. Since the position of the electrometer relative to the 
circumstances was not changed during the observation, the observed value is regarded 
as the relative value of the electric field at the point O. . The instrument used for 
mesuring the space charge is a cubic brass cage, whose side is 90 centimeters long, 
and at the center of which a radioactive collector was held by a bamboo rod. The 


scale value of this instrument is, 
0=0.58 Ve.s.u./m? 


where V is the electric potential of the collector measured in volt. 

Among the results observed, the most reliable values at point O are inserted in 
Fig. 2. On seeing this, it is evident first that the electric field under the smoke is 
always positive and also the space charge is almost always positive except in the short 
durations in which it has a small negative vaiue. Secondly, the variation of the 
electric field and the space charge are parallel. For example, dividing the data into 
the groups A, B, C, D and E, as shown in Fig. 2, we obtain the correlation coefficients 
between the electric field and the space charge in each group, which are, 


Group Sample number _ Correlation coefficient 
A 39 0.56 
B. 70 0.46 
G5 74 0.63 
_D 35 0.22 
E 98 0.39 


In groups A, B and C, a good positive correlation is observed, but not so in groups 
D and E. However, the electric field strength and the space charge have a similar 
a in small variation, even in groups D and E, only a linear relation between 
them is partly broken, owing to the superposition of large variations, which are not 
related linearly to each other. We notice thirdly that the positive maximum values 
of the space charge fluctuate pretty largely but its negative minimum values remain 
fairly constant, that is, near 0.5 es.u./m’. Further, it is remarkable that the space 
charge was generally in inverse proportion to the distance from the center of the 
smoke, and that the distance where the space charge changes from positive to nega- 
tive, is found to be about 60 meters, as shown in Fig. 7. The distance between the 
center of the smoke and the observing point is estimated with eye-observation. In 
practice, the horizontal and vertical distances "9.€ Z) were observed with respect to 
a co-ordinate, having an origin at the observing point and the X-Z plane being at 


12 
right angles to the center line of the smoke cylinder. The eye-observation of the 


distance may have an error of +10 meters. 


3. Charge Distribution In and Around the Smoke 

On the basis of the observed data stated above, we consider the charge distri- 
bution in the smoke. It must be. noted first that, since the electric field is always 
positive at the ground surface, the main part of the smoke may be expected to have 
positive charge. The intensity of the field strength is estimated ten times as large as 
the ordinary value, considering the sensitivity of the electrometer and its relative 
position to the surroundings. The observed value of the space charge frequently 
became 3 e.s.u./m’. The smoke remains in general the shape of circular cylinder 
without being broken to pieces. Taking these facts into consideration we may assume 


that the charge density of the smoke p) has the following expression, 
a2 


we rae 3 
pt) =0n( @ ag ee r<a 


Sea) a Y>a 


(1) 


where 7 is the distance from the center of the smoke cylinder, and pm, ¢ and a are 
constants which are estimated respectively, fed i 
(m=) e.s.u./m§, o=25m, ~a=70m. 

- The space ‘charge p, which is actually observed, is the sum of 00 and the charge 
induced by the smoke. Since the mobility of the smoke particle is considered as 
smaller than that of a large ion, 0) may be assumed to be constant within a compara. 
tively short time after the smoke starts from the crater. The distance between the | 4 
crater and observing point is about 350 meters and the wind velocity was about | 3 
meter/sec. during the observation, so the smoke takes: about two minutes on the — 
average to ESS yer the distance fre 
50 Imin.t0 20 30 6080 2ninl0 030 40 50 Iain 0 20 30 erving pc "4 

3 inducing negative charg g al 
a BESS aig * ped de 3 


ne The ‘schematic aspect of the’ inducing 
_ negative —— prommte stnceitokes 


— 
w 


becomes about 0.06 coulombs and is She 


or two orders less than that expected cal ze 

on Volcano Asama and Volcano Yake. 3 “% 

It may be natural, since the smoke of 2 § es 

Volcano Azuma is fairly weaker than Sv. ou 
2 


that t untains. 
hat of the above two mountains. Now BNO oe! dae seri me or nc 73 


assume that the positive and negative Tim) 

Fig. 4: The space charge pp and the electric 
field intensity Ey of the smoke in 

nitial state outside the smoke, then the the initial state. 


ions are spread homogeneously in the 


positive ions move in the direction of 7. Let us consider first the motion of a positive 

_ charge outside the smoke. Using the equation (3), we have, ; 

dr*/dt=kE,=2kCQ/r*, rt a (4) 

where Q is measured in practical unit and C is a constant, having a value of 9x10", 
_ Integrating the above pacer we obtain, 


‘rts (4kCQt+s) (5) 
"where s is gifs position of the ion considered at the initial state. se ag this, we get, 
dr*/ds = S/rt, (ay 
So that, : 
— Co*(h) dk) =Ce f (R) dR je=e =Cot(k) dk )e=0 =const. : (7) 
ris valid in ‘the region, in which the positive charge of mobility k can exist. abe 
zs eee the charge of an ion f “<r e 
: - (dk : spectral density function of positive ion, whose mobility is between k 
fis i. a and k+ dk. ~ ec tse ata aga aie merectce d 


ns paently, the charge density with soni po positive ions at a certain place and time 


f >: 


P a rd 

ae rhe =f c am! i : aT } 7 a P = = ’ B 

Ee Ses Pye Fi kydke ot aos (8) 

935 Iuserrrot a etd On “Cee 

2ae a ‘ = ma eee : Sky 

Wate ore one another’ eee Bc: 39 diet eobeastas vee he 
hi Fis elie j kas st rf K. zat ~ = Ea ACQt Vetis Ppp tan 5 Ess efi 2 dds 3 eS; o 


14 
where, X=log he 
and i=1: small ion 


2: intermediate ion 
3: large ion. 
For simplicity, let us consider the case i=1, 


In this case f (R) dk becomes, 
. Ff (k) dk=216 exp (—€) dé 
&=1.543 log k 


Making use of the above density function we obtain, 


)=—307x10-" {4 ee ae | 
0 (%, Sing : tele ee 


Tw JO 


where, 


£3 log | eit oc a 
aie: onl 


A, =334, /1=0,... - ¥y=0,458. 


Next, we must consider the case. r<a. We found easily, however, that the charge dis- 
tribution of the negative ion inside the smoke i = identical with. that outside the imal in 


die 


= aes 3 sec, the first approximation, aft- 


lar to the case of outside 2 the 


‘ phifey 
smoke. “Thus: we have, 
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le . dit * < , x 
place. I’o consider the cases where i=2 or 3 will also serve to enlarge the density 
of induced negative charge. Next we ‘consider the electric field under the smoke. 


Taking the image charge against the earth’s surface into consideration, the electric 


field on the ground’ surface is given by, 


Ex=0, 2=—(4ZCQ)/R? 
Volt/cm 


where, 


YR? = Dee Se. (16) 


is the distance between 


the center of the smoke 


and the observing point. 


Using meter for the unit 


of length, we can rewrite 


the above equations as, 


E2= — 2120 2/R? 
Volt/em (17) 


The electric field caiculat- 


ed with the above equa- 


tion is shown in Fig. 6. 


O- au 200 0 50 60 70 80 9 BOX Gn) 


Fig. 6: The contour line showing the position of the center of 

the smoke which gives the same electric field intensity 
how::the charge distribu- ,at the origin (The value on the contour shows the elec- 
tric field intensity at the origin measured in volt/cm). 


4. Discussion 


We shall examine 


tion obtained in the pre- 
vious section will explain the observed results. From Fig. 5, it is expected that, as 
the center of the smoke moves, the maximum value of the space charge considerably 
fluctuates in the positive region, but in the case of negative region, its minimum 
value may have an almost constant value. Since Fig. 5 is that calculated only for the 
stall ion, the effects of the intermediate and the large ion must be taken into account. 
The calculated value of = (7) in the case that the effect of the intermediate ion is 
considered, is illustrated in Fig. 8. The depth of the: channel of negative charge 
around the smoke is determined by the shape of the ion spectrum. It is likely in the 
actual case that the steplike part of the charge distribution in Fig. 3 will be made 
more smooth by the turbulence of the wind. In either case, the depth of the channel 
of negative charge around the smoke may be nearly constant within a- comparatively 
short space of time. The problem we must consider next is whether we are justified 
in neglecting the effect of the ground in equation (4). We may approximate it more. 
closely, using the equation (17) instead of (3). Speaking qualitatively, to make use 
‘of the. equation (17). results in broadening the width of the channel of negative charge 
just under the smoke twice as large as that just above it. In the actual case, how* 


ever, po w” will be distorted considerably as the smoke approaches to the ground, so 


16 
that detailed calculation with the equation (17) may be rather insignificant. 

As a transient case we introduced the equations (14) and (15), that is, we en- 
tirely neglected the recombination of the ions. But, if we consider it, the result 
obtained may not be modified so largely in this case because the time after the smoke 


starts from the crater is comparatively short. 
Though it is such a rough approximation, the charge distribution obtained 


X,2.R 
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Fig. 7.: The comparison of the calculated values of space charge and the potential gradient 
with those of the observed values. 


above can explain the various phe- 


2 
nomena shown in the second section: Se 

From the equations (15), (16) and (177 § . 

it is clear that the variation of the S —7(m) 
space charge and that of the électric field x 0| YB. 

~ 
strength are almost parallel to each & \ 
Ss small ion 


y 
= 


other, both inversely proportionate to 
R. The values of o and E obtained from 
actual observation and those by calcula- 


? Ui , tf. t o 
Fig. 8; The curve of p(7) in the case that 
both small and intermediate ion are 


tion in which the case ¢=140 sec. is considered. 


adopted, are illustrated in Fig. 7. So far as the smoke does tot approach very near 


to the observing point, the calculated values sufficiently coincide with the observed: 
values. Further, the smoke often changed its density, the state of which is roughly 
shown in Fig. 7, and the wind velocity was not always constant, therefore the width 
of the channel of negative charge is naturally expected to fluctuate, but a detailed. 
discussion of such effects is beyond our present treatment. 

We should like to express our hearty thanks to the members of ae of RS 


7 


Fukushima Meteorological Station for their great help in this observation and also 

to Dr. H. Hatakeyama, Director of the:Meteorological Research Institute, for his kind 

suggestion and encouragement in the present research. a1 
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The Annual- and Diurnal Variations of Cosmic-Ray Intentsity 
and the Temperature Effect 


By Kaichi MAEDA and Tomoshige SUDA 


Meteorological Research Institute, Tokyo 
Abstract 


In order to estimate the influence of the diurnal variation of 
the atmospheric temperature on that of cosmic ray intensity, we 
investigated the bihourly data of cosmic ray intensity at Huancayo 
and aerological data at San Juan. As the mean atmospheric tempe- . 
rature, the altitude differences between 100 mb. and 1000 mb. levels 
were used. We ascertained that the amplitude of the diurnal change 
of cosmic-ray intensity at Huancayo exceeded 0.5% at least (i.e. twice 
that of the observed value), if the temperature effect was taken into 
accounts. 

1. Introduction. 

Previously, we disccussed in brief (1) the influence of the diurnal variation of 
mean atmospheric temperature upon that of cosmic ray. The upper-air data used 
in that discussion, however, were not accurate enough to treat of such a problem _ 
quantitatively. Therefore, we have investigated on this problem with respect to the J 

- cosmic ray data of Huancayo, taking reliable aerological data of the diurnal change 


Apt Jun. a Oct. i into accounts. B this paper, the tem: 
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high latitudes than in low, and this is of course due to the so-called. temperature effect 
ont he cosmic ray. Therefore, in order to. estimate the coefficients of this effect for 
these data, the annual variations of the altitude difference between 100mb. and 
1000 mb. levels at Swan Island (15° N), San Juan (17° N), Washington (38.°9 N) and Thule 
(77 °N) are evaluated as in Fig. 2. The curves are computed from the daily upper 
air data published by Weather Bureau, 
U.S.A. in 1950. 


As a rough estimation of the 


Km 
164 


16.6 
coefficients, taking the ratio of both 


os these amplitudes, it follows 


15,2 -9.0%/km for Huancayo 


55%/km for Cheltenham 
and -3.4%/km for Godhayn. 
Thule (77'n ) Among these values, -5.5%/km is 


15.6 


only plausible because Washington and 


160 Cheltenham are very nearly situated 


and no complete upper air soundings 
Fig. 2 The Annual Variation of the Altitude 
Difference between 100mb. and 1000mb. 
Levels. Washington. Moreover, the value -5.5 


for the year are available except for 


%/km is exactly the same with that determined by the calculation of regression 
coefficient. Dashed lines in Fig. 1 denote the curves corrected by these coefficients 
using the curves in Fig. 2. In the corrected curves, the semi-annual changes are 
predominant, but this is merely apparent because both data are taken from different 
years. Zi eae 

3. The Diurnal Variation of Cosmic-Ray Corrected for the Temperature Effect. 

In order to consider the temperature correction which affects the forms of the 
diurnal variation of cosmic ray intensity at Huancayo, the diurnal variation of the 
altitude difference between 100 mb. and 
680 mb. levels is calculated as the low- 
est curves in Fig. 3 from the data of 
three-hour interval observations by 
radio-sonde performed at San Juan 
during October and November in 1944, 
for the data of Huancayo are corrected Atstude Difference 

; : : (Joo-680 mb) ~. - 
to the values of 510 mmHg (i.e. 680 mb) 
for the pressure effect. Usually the 
day-time temperature in upper atmos- itcne 


phere is overstimated owing to the Fig. 3 -The Diurnal Variation of the Cos- > 
; me : if _ mic-Ray (upper) and the Altitude. 
solar heating of meassuring instrume Differnce between 100mb ae 680mb 


and this gives an enlarged amplitude Levels (below). 
of the diurnal changes of the isobar heights in the upper atmosphere. But, in. aiid 
‘respect, those results are corrected cautiously by the graphical method (3) and per: 
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mit us to use them for our purpose because the diunal variations of aerological states 
may be regarded to be similar in San Juan and Huancayo though they are ina diff- 
ence hemisphere from each other. E ? Yj 

Now, to estimate the temperature effect for the date of Huancayo, let us consider 
the altitude dependency of cosmic ray intensity in the lower atmosphere simply as — 

I=, exp{—B (h—z)} 
where I is the intensity at the height zkm above the sea level, J, is the intensity at 
the meson producing layer whose height is assumed to be 4km above the sea level, 
and # is the mean decay coefficient of the meson component and this is equal to the 
reciprocal of L, the mean life range, wlich is expressed by 
Le, | 
mp 
if we indicate by Tt, the life time of meson (nearly 2y sec.), 
by m, the rest mass of meson (nearly 10% eV/c’) 

and by ~, the momentum of meson in MeV/c. 
Then the coefficient of temperature effect at sea level, b, is written by 


fale Sat B Cola 
b= an 1/L 
; 1 OL 
. (fea Bl I 
and at the height of z b= T oh LL 
2=Z “y 
Thus, we get eal E PE aS 7 Yet 


yes L aad L’ represent the mean life range at 2-0 and 2=Z respectively and. ace 
cording to Bennedetto (4) L=14.2 km (this corresponds to p=2.7x 108 MeV/c). Iti is. 
clear that L’ is shorter than L, for the momentum Spectrum increases with the ale 
3 titude more rapidly in lower regions than in higher. Put L’= 8.7 km (this cor res) onds 
— to p= 1.7x103 MeV/c), then we get J’ =—9.0%/km as ahowe above, but thi is v 
seems to be too large, so that as the lowest estimation we take 5.5%. km. 5 


é rected curves by this factor is shown in 1 Fig. 3-by the dotted Heeb, the 
indicates the uncorrected variation. > ee a eee ee aries 
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Our thanks are due to Dr. K. Takahashi, Forecast Research Laboratory of this 
Institute, for his kind informations about aerological data. 
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LETTER TO EDITOR 


A Newly Designed Differential Electrometer and its 
Application to the Simultaneous Measurment of Air 
Earth Current and Potential Gradient 


F.J. Scrase made a simultaneous measurement of air earth current and atmos- 
pheric potential gradient using radioactive collector and Wilson’s test plate. In order 
to eliminate the induced charge on the test plate due to atmospheric electric field 
and to get pure air earth current, Scrase employed the quadrant electrometer differ- 
entially. In this case the electric potentials on each two quadrants are same sign, 
so that by suitable adjustment two potentials are compensated each other, and the 
electrometer is being charged by pure air earth current only. When the sign of two 
potentials, which are to be compared each other, are opposite, we can not use quad- 
rant electrometer differentially. For this purpose new type of differential electrometer 
becomes neccessary. This new electrometer is designed and constructed by the author. 

The principal parts of the electrometer consist of, as shown in Fig. I a, two 


needles fixed perpendicularly each other to the metal stick, on which a small mirror 


_ MECHANICAL 
_. COLLECTOR ‘- 


POSITIVE TEST PLATE 
FIELD ee 
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is attached at the middle point, and the double quadrant system are located above 
and below. As an example, the author applies this electrometer to the simultaneous 
measurements of air earth current and potential gradient using rotating collector 
and test plate. In the case of this collector, rotating vane is earthed when exposed 
to the atmospheric electric field, and connected to the electrometer when covered 


by the earthed shield, so that the sign of the potential obtained is negative when the 


atmospheric electric field is positive, or vis versa. 
Fig. I b shows the connection circuit. The electrometer E, should be deflected 


_ by the charge accumulated by the air earth current flows into the test plate. The 


3 


charge induced on the test plate by the atmospheric electric field is compensated by 
the action of rotating collector through variable condenser adjusted suitably. In this 


_ figure the manner of compensation is shown in the case of atmospheric field is posi- 


a ae he wpa ver 


; 


: 


Ki 4 is Pee 
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tive. We can find the fine structure between air earth current and potential gradient 
by this method. Employing such a electrometer some measurements were carried 
out. The detailed results of this measurements will be reported in the near future. 
The author expresses his thanks to Dr. Prof. M. Hasegawa and Dr. Y. Tamura 
for their valuable advices. 
| By Miyoji Goto 
Aso Magnetic Observatory, Kyoto University 


(1) FJ. Scrase, Geophys. Mem.,No. 58, 1933, 
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